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We present a detailed study of pion and rho mass, decay constants and quark mass in Lattice QCD with two 
flavors of dynamical quarks. We use Wilson gauge and fermion action at j3 = 5.6 on 16 3 x 32 lattice at eight 
values of the Wilson hopping parameter in the range 0.156 - 0.158. We perform a detailed investigation of 
the effect of gaussian smearing on both source and sink. We determine the optimum smearing parameter for 
various correlators for each value of the Wilson hopping parameter. The effects of smearing on observables are 
compared with those measured using local operators. We also investigate systematic effects in the extraction 
of masses and decay constants using different types of correlation functions for pion observables. We make 
interesting observations regarding chiral extrapolations and finite volume effects of different operators. 
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I. INTRODUCTION 

Among the various formulations of Lattice QCD, the Wilson formulation [ 1 ] (both gauge and fermion action) proposed long 
ago is still attractive due to many reasons, first and foremost being the conceptual simplicity. The Wilson formulation preserves 
discrete symmetries of the continuum formulation which simplifies the construction of lattice operators that correspond to the 
observables in the continuum theory. However, because of the explicit violation of chiral symmetry by a dimension five kinetic 
operator, Wilson formulation has been difficult to simulate at light quark masses. Absence of chiral symmetry implies that the 
"physical" quark mass is no longer proportional to the bare quark mass and the quark mass renormalization is no longer only 
multiplicative. Lack of chiral symmetry means that the Wilson-Dirac operator which is a sum of anti-Hermitian and Hermitian 
terms is not protected from arbitrarily small eigenvalues and may lead to zero or near zero modes for individual configurations. 
This is the infamous problem of "exceptional configurations". This leads to convergence difficulties for fermion matrix inversion 
which is an integral part of both the generation of gauge field configurations with dynamical fermions and for the computation 
of hadronic correlation functions via quark propagators. This poses difficulties for lattice simulations with Wilson fermions in 
the chiral region. The situation has improved recently partly due to the finding [2 1 that the numerical simulations are safe from 
accidental zero modes for large volumes. 

Very generally, the major challenges of probing the chiral regime of lattice QCD are: a) to be able to achieve small lattice 
spacing a to reduce scaling violations, b) to have small current quark masses to have reliable chiral extrapolations, and c) to 
have large enough physical volume of the lattice to avoid finite size effects. In addition, lattice QCD investigations suffer from 
uncertainties regarding the determination of the lattice scale and inaccuracies in derived quantities of interest, e.g., hadronic 
masses, decay constants etc. Many of the above issues are connected. 

We have planned on a detailed lattice QCD investigation with Wilson and Wilson-type fermions. In this work 
we have taken the first step towards addressing the above issues systematically. We study smearing of mesonic op- 
erators with standard (unimproved) Wilson gauge and fermion actions with 2 fully dynamical light quark flavors on 
16 3 32 lattice, j3 = 6/g 2 = 5.6 at 8 values of bare quark masses (corresponding to fermionic hopping parameter K — 
0.156, 0.1565, 0.15675, 0.157, 0.15725, 0.1575, 0.15775, 0.158). The lattice scale reached is respectable and a ~ 0.08 
fm (or 1 ~ 2.45GeV). It is important to have a small enough a to hopefully have small scaling violations and then evaluate the 
low lying spectrum of QCD accurately and to study cleanly chiral extrapolations in finite volumes. We plan on having larger 
volumes and larger /3 in the near future with more sophisticated algorithmic developments discussed in [2 3[ffl[5]]. In the present 
work our emphasis is to have fully dynamical simulations at an extensive set of K values and study accurate determinations of 
the masses and decay constants in as many ways as possible. 

In a separate paper (6), we investigate quite elaborately the determination of the lattice scale from the potential between a 
heavy quark-antiquark pair. We study the observed change of ro /a with bare quark masses where ro is the Sommer parameter 
Q and resulting uncertainties of scale determinations and also discuss its relation with the chiral extrapolation of observables. 
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For a recent review on various approaches to scale determination and associated issues, see Ref. [8]. 

The masses in Lattice Gauge Theory are calculated from the asymptotic behaviour of Euclidean time correlation functions. 
The contribution from the lowest mass state dominates for large time. Thus, to get a clean signal one has to make the length in 
the time direction as large as possible. But this is computationally expensive. Furthermore, as the time increases, the signal to 
noise ratio gets smaller which results in larger statistical errors. To reduce the statistical error, we need to increase the number of 
measurements (configurations). This is also computationally expensive. Hence there is a practical need for techniques that allow 
one to reliably extract observables still working on moderately sized time direction and not too large number of configurations. 

For the measurement to be in the scaling region, one needs to go to smaller and smaller lattice spacings. But as the lattice 
spacing decreases, the physical hadron state will extend over more and more lattice spacings. Thus a state created by a local 
operator from the vacuum will have less and less overlap with the physical state as the lattice spacing decreases. To improve 
the measurements, it is useful to use operators that have larger overlaps with the physical state (which is extended). This will 
naturally lessen the contamination from higher mass states to the correlation functions. Smearing of operators is one way to 
achieve this goal, thereby reducing the necessity to have both large time direction and large number of configurations. 

One of the earliest references discussing the need for smearing is Ref. [9|. Earliest suggestions for smearing involved Dirac 
delta functions on source time slices, the so-called wall source. See for example, Ref. ifTUl . So far, in the literature, smearing has 
been implemented in two different methods, namely, the gauge invariant method and the gauge fixed method. Gauge fixed wall 
source smearing, exponential smearing and gaussian smearing belong to the second method. Gauge fixed wall source smearing 
is used in some work of JLQCD ifTTI and MILC 1121 Collaborations. Exponential smearing is used in some works of CP-PACS 
[13 1 and JLQCD [14| Collaborations. Gaussian smearing was first considered by DeGrand and Loft 1151 . and used, e.g., in Ref. 
USE)) (with dynamical quarks) and in Ref. Ifl8l [191 l20l |2T1 (with quenched quarks) in different contexts. A clear discussion 
of smearing (in the context of gaussian smearing) and the associated Fourier Transformation "trick" is provided in the Ph. D. 
thesis of Hauswirth |[T6l . 

To the best of our knowledge, a systematic study of the effects of gaussian smearing on source or/and sink operators does not 
seem to be available in the literature. In this work we carry out a detailed investigation of the effect of gaussian smearing on the 
source or/and sink meson operators. 

At a particular value of the smearing size, a variety of correlation functions are to be measured to extract meson observables. 
For example, for the pion, PP, AA, AP and PA correlators can be used since both P and A carry quantum numbers of pion. (Here 
P and A denotes psudoscalar and the fourth component of the axial vector densities respectively.) Mass gaps and amplitudes 
(coefficients) are extracted for the smearing radius .so ranging from 1 to 8 in steps of 1. We also carry out extensive calculations 
with local operators. This helps us to quantify the systematic effects in the determination of masses and decay constants. 

We present results for the pion mass, the rho mass, their decay constants and the quark mass (all in units of the lattice constant 
a). Although this paper is not intended to deal with the issues of chiral extrapolation of pseudoscalar observables, especially 
at our relatively large masses with 16 3 32 lattices, we make interesting observations regarding consistency with %PT. In this 
connection, we also comment on finite volume effects on part of our data at the largest K. We have also noticed varying finite 
size effects on different pion operators. 

This paper is organized as follows. In Sec. [Iljwe present the lattice action and the details of the simulation. Expressions for the 
local observables that we calculate are given in Sec. [Ill] Sec. [jVjintroduces the gaussian smearing and different fitting ansaetze 
for the correlator data analysis are given in Sec. [V] Details of the implementation of smearing are in Sec. VI Results for pion 



and rho observables are presented in Sees. |VII| and [VIII respectively. Finally, the summary and the conclusions are presented 
in Sec. [IX] For the sake of completeness and clarity, the Fourier transform method employed in the case of sink smearing is 
presented in Appendix |A[ 

II. SIMULATION 
A. Action 

We have performed simulations with the standard Wilson action S = Sf + Sg where the standard Wilson gauge action is 
given by 



1 

1 — - Re Trace Up 



with /3 = (g is the 51/(3) gauge coupling), the elementary plaquette Up being the product of SU(3) link fields U around the 



elementary square of the hypercubic lattice and the standard Wilson fermion action is given by 

S F [y,\jf,U] =Y,Wx M xyYy with (2.1) 
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FIG. 1 : Average value of the Plaquette versus K. 



M xy = 8 xy — K (r-Yn) Ux,n $ 



'x+ti,y 



(2.2) 



We have suppressed spin, color and flavor indices in the quark fields yf and yf and the fermion matrix M. As usual, we have 
taken the Wilson parameter r = 1. We consider two degenerate light quark flavors, i.e., Nf = 2. 



B. Details of Simulation 

The gauge coupling j3 = 5.6 and the lattice volume is 16 3 x 32. The hopping parameter K"=0.156, 0.1565, 0.15675, 0.157, 
0.15725, 0.1575, 0.15775 and 0.158. At each K we have generated 5000 equilibrated configurations with the standard HMC 
algorithm (with even-odd pre-conditioned Conjugate Gradient for inversion of M^M) and performed the mesonic correlator 
measurements separated by 25 trajectories. The time step size is chosen to be 0.01 and the number of steps per trajectory is 
100. The stopping criterion for the residue of the conjugate gradient y] R \ 2 / \ % | 2 < 10~ 5 is used for all K values except 0.158 
(for which it is 10~ 6 ) where R = M^M^ — % (% and £ are respectively the source and the iterative solution). Stabilized BICG 
algorithm is used for the inversion of M needed for the mesonic correlators, the stopping criterion of the residue for very accurate 
inversion in this case is chosen as 10~ 9 and 10~ 10 . 

The parameters are chosen such that the acceptance rate is above 75 % for K less than 0.157 and around 65 % for K above and 
including 0.157. 

Fig. [T] shows that the expectation value of the average plaquette has a smooth behavior in dependence of K. 



C. Auto correlation 



For any particular observable 6, autocorrelations among the generated configurations are generally determined by the inte- 
grated autocorrelation time for that observable. For this purpose, at first, one needs to calculate the unnormalized autocorre- 
lation function of the observable measured on a sequence of equilibrated configurations as 

C ff {t) = -i-^Ww K + , - (<?)«) (2.3) 
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TABLE I: Integrated Autocorrelation time T, )M for Average Niter and Plaquatte as a function of K. 



where, instead of the usual single mean-value estimator, we have used the "left" and "right" mean-value estimators 



1 



N-t 



L = 



N- 



l 



N-t 



R = 



,v=l 



N- 



7 5 Gs+t 



(2.4) 



for a faster convergence as claimed in 1221 . 

Following the "windowing" method as recommended by Ref. 11231 . we have calculated the integrated autocorrelation time as 
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(2.5) 



where 



is the normalized autocorrelati on fu nction, 
l 



p ff (t) = C^(t)/C e 



(0) 



(2.6) 



the 



The factor of i in the Eq. (2.5 1 is purely a matter of convention and is chosen such as to guarantee that for N - 
normalized autocorrelation function behaves as exp (—t/T? t ) with 3> 1. 

If plotted against the variable cutoff f cut , the resulting values of T^ t ideally show a plateau but instead a peak or monotonous rise 
is also observed sometimes. Following the suggestions of [23], we have drawn two straight lines f cu t/4 and f C ut/10 respectively 
and have chosen inspecting only the segment of the curve lying within its intersections with the two straight lines. This 
ensures a balance between the noise and the bias in the estimation of the autocorrelation time. 

We have measured the integrated autocorrelation times for the average values of the Plaquette and the number of Conjugate 
Gradient iterations, called Niter, during HMC trajectories for the generated configurations for each value of K. The errors are 
calculated by the single omission jackknife method. Our results are presented in the Table|I] Despite fluctuations, T; nt in general 
increases with larger K especially as K approaches K c . 

A typical example of the expected plateau of T; nt as a function of f cut for both the average Niter and Plaquette at K — 0.157 is 
shown in Fig. [2] 



III. OBSERVABLES WITH LOCAL OPERATORS 



We measure the charged pion and the rho propagators to extract masses, decay constants and the quark mass at each K. 
For pion, we measure the following zero-spatial-momentun correlation functions on a L?T lattice as functions of the Euclidean 
time t : 
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FIG. 2: Integrated autocorrelation time v ml for both average Niter and Plaquette at K=0. 157. 



00 = PP or AA and 0\02— AP or PA where P = qiYsqj and A4 = ^,7475^; denote the pseudoscalar density and fourth 
component of the axial vector current (z and j stand for flavor indices for the u and d quarks, for the charged pion i ^ j). 

For clarity and to set up our notation, in the following we discuss all the possible ways of determining the decay constants 
and the PC AC quark masses. 



A. Pion Decay Constant 



The pion decay constant F n and the quark mass m q from PCAC or the axial Ward identity are respectively defined, in the 
continuum, via 



(0|A M (0)K/>)) = V2F n p^ 
dfiA^(x) = 2m q P(x). 



(3.1) 

(3.2) 



Since we measure the PP, PA, AP and AA correlators, we have a variety of ways to compute the pion decay constant and the 
PCAC quark mass. 



Method I: From the AA correlator 



C(t) = (0\Al(t)A 4 (0) |0) 
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(3.3) 
(3.4) 



where C M = ^- | (0 | A 4 (0) | n) | 2 , and the pion decay constant as defined in Eq. (13 .1 



follows 
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The factor 2k above accounts for the difference in normalization between the continuum and the Wilson lattice fermion actions. 

It is obvious that the m K that appears in Eq. ( |3.5[ ) is numerically obtained from the AA correlator. If numerically m K is obtained 
from another correlator, say the PP, then one needs to put that mass in the fit Ansatz for the AA correlator to obtain the coefficient 
C M . 

Method II: From the PP and the AP propagators 
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Similarly using the PP and the PA propagators 
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(3.6) 
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For numerical evaluation of F K using Eqs. (3.8 i and (3.9 1, the values of m n computed from PP, AP and PA correlators have 



to be the same, which is hard to achieve numerically. Hence it is advisable to use the best determined mass from one particular 
correlator and then use that value in the other correlators to determine all the coefficients C. The same applies also to the 
determination of the different PCAC quark masses discussed below. 



B. PCAC Quark Mass 

Method I: From d^A^ (x) = 2m q P(x) , summing over spatial coordinates, 

£^ A ftM = 2m q 1 £P(x) d A A A {t) = 2m q P{t). 



(3.10) 



Taking the matrix element between the vaccum and physical pion states we have 

m,(0|A 4 (0)|w> = 2jfi,<0|i'(0)|ff) 

which leads to 
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(3.11) 



(3.12) 



Method II: Using PCAC 



Summing over spatial coordinates 
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(3.13) 
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At large t , 



which leads to 
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C. Rho Mass and Decay Constant 

There are two different definitions of p deacy constant used in the literature. 

m 

<0|V M (0)|p) = e M (3.17) 
Jp 

where f p is dimensionless and 

(0|7 M (0) |p) = e^V2F p m p (3.18) 

with V^(0) = ^(0)7^^(0) and is the polarization vector of rho. 
Here F p is dimensionful and F p /m p = (fp) . 

Mass and decay constant of rho are calculated from the correlation function 

C(t) = (0 1^(^3(0)10) (3.19) 
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where, C vv = ^ | (0 | V 3 (Q) | p) | 2 . Thus 



c vv 

Fp/mp = l/f p = 2k J— . (3.21) 



We compute the decay constants and the quark masses according to the expressions given above, but to obtain their values in 
the continuum one needs to multiply with appropriate factors of renormalization constants Zp, Z& andZy associated with lattice 
pseudoscalar, axial vector and vector densities respectively. Although we have made approximate estimations of Za and Zy in 



Sec. VIII we do not quote numbers in the continuum in this paper. 



IV. GAUSSIAN SMEARING 

To increase the overlap with the hadronic ground state, many, if not most, QCD spectrum calculations traditionally use the 
method of smearing the hadronic interpolating operator, essentially making the hadronic operator spread around their central 
location in space. In this work, for the pion and the rho operators, we have used the so-called gaussian smearing where one uses 

a shell model trial wave function with one variational parameter, <p(r) ~ exp ^— (r/so) 2 ^j . For details please see Ref. f |[T5l l. so 

is the smearing size parameter. 

An advantage of the gaussian smearing is that the smearing function, for example for a meson operator, separates into two 
factors one belonging to the quark and the other to the antiquark. This has certain numerical advantages, especially for sink 
smearing. However, the smeared operators are no longer gauge-invariant because the quark and the antiquark are spatially 
separated. We have used Coulomb gauge fixing to obtain non-zero expectation values. 



The coefficients C ff[ 61 for local operators are needed, as described in Sec. Ill to calculate the decay constants and the PCAC 
quark masses. One can take a mixed approach where one evaluates the hadron mass from a hadronic propagator which uses 
smearing at either the source or the sink or at both places, and at the same time calculate the coefficient C ffi 62 from the local-local 
propagator (as pursued in Ref. 02)). For this to work, large time-extents are necessary. However, if one wants to calculate the 
local-local coefficients from the smeared propagators, one needs to calculate the propagators with all combinations of smearing: 
i) local sink and smeared source (Is), ii) smeared sink and local source (si), and (iii) smeared sink and smeared source (ss). If 
they produce the same hadronic mass at large euclidean times, by combining the coeffecients of the three, all calculated with the 
same smearing parameter sq, one is able to calculate the coefficient corresponding to the local-local propagator. This is done as 
follows. 

The large Euclidean time behavior of the correlation function involving local (unsmeared) operators is jz/ and 38 

(0 I < inkt (f) ^ source (0) I 0) ^ Cf( m e- m »> 
where, Cf( m = j^(0 | ^ sinkt (0) \H)(H\ S§~(Q) \ 0) . 
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Assuming the lowest mass obtained at large Euclidean time to be the same for each of the Is, si and ss correlators, it follows 
that 



Qi — Ci s C s i/C ss 



(4.1) 



V. ANALYSIS OF MESON CORRELATION FUNCTIONS 

For our initial investigations we performed four types of fits using exponential functions for the correlation functions of the 
pion and the rho. For example, for the pion propagator, the fit ansaetze involving two, three, four and five parameters are as 
follows. 
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(5.1) 

In the presence of sea quarks, pair creation from vacuum becomes possible and the creation of two quark-antiquark pairs 
results in two extra pseudoscalar mesons. Thus in the presence of sea quarks, the next higher state can be expected to be a three 
pion state [3 1. The second and the fourth fit ansaetze above are motivated by this physical picture. The rho propagator can also 
be treated in the same way, the next higher state there being 2m K + m p . 



VI. DETAILS OF THE IMPLEMENTATION OF THE GAUSSIAN SMEARING 



As mentioned earlier, in this paper we present results with two flavors of fully dynamical Wilson (unimproved) quarks (with 
standard plaquette Wilson gauge action) at j3 = 5.6 on 16 3 32 lattices at an extensive set of the fermion hopping parameter, viz., 
k = 0.156, 0.1565, 0.15675, 0.157, 0.15725, 0.1575, 0.15775, 0.158. 

At this j3 and lattice volume, similar Nf = 2 calculations have been done before at a few K values (for example, see Ref. 
Il22l ). The previous results help us to have a cross-check and have belief in our numerical procedure. Some data are available 
at larger volumes with the rest of the parameters staying the same; these give some indication of the finite size effects in 
our results. However, we like to mention that, to the best of our knowledge, ours are the first calculations at j3 = 5.6 at 
K = 0.15675, 0.15725 and 0.15775 for any lattice volume. 

At each K, we have investigated in detail the PP, AA, AP and PA correlators for the pion and VV correlators for the rho. For 
each operator, we have used source smearing (Is), sink smearing (si) and smearing at both source and sink (ss) (see Appendix [A] 
for the Fourier transform method for the smeared sources). Then at each K, with each operator and each sink-source smearing 
combination, we have tried eight values for the gaussian smearing size parameter, viz., so = 1 to L/2 — 8 (in increment of unity) 
in order to have optimum smearing (please note that so —> 00 produces the so-called wall source). 

For comparison, we have also investigated all the above correlators without any smearing, i.e., using local sink and local 
source (denoted with subscript 11). 

With or without smearing, our general strategy was to determine the pion mass from the PP and the AA propagators only 
and not from the asymmetric propagators AP or PA, by reaching a plateau with respect to f m ; n (minimum t used for fitting the 



formulae 5.1 In general, the quality of the mass-plateaux for the AP and the PA correlator was not of the same high quality as for 
the PP and the AA correlators (AA noisier than PP), more true with smearing. As a note of caution, let us also mention that at a 
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FIG. 3: Pion mass obtained from unsmeared PP correlator at K" = 0. 157 with four different fits as function of /„ 
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FIG. 4: Pion mass from one exponential fit (two parameters) as a function of sq. 



few large K values the local-local AA correlators were very noisy, resulting in inaccurate and perhaps unreliable determinations 



of the all relevant quantities as apparent from fourth data-column of Table III 

We have used the well-known technique of effective masses by taking ratios of propagators for adjacent time-slices, but only 
as a rough guide (especially for smeared propagators). We shall not show any effective mass plots in this paper. 

In Fig. [3| w e show, for the unsmeared (11) PP correlator, the results at K = 0.157 for m K obtained with the four different fit 
Ansaetze |5.1| as functions of f m j n (here and in the subsequent figures, tmi n /a in the abscissae is indicated as just f m ; n . Firstly, 
we note in Fig. [3]that all the different fits produce good plateaux for appropriately large f m ; n and they all agree with each other. 
The plot also shows that 3m K is the next higher state. Unless otherwise specified, all errors in this plot and others to follow are 
single-omission jackknife statistical errors calculated from 200 jackknife bins of correlator data. 

With unsmeared propagators, although it is obvious that longer time-extent would have helped, we have been quite successful 
in obtaining a reasonable plateau even at lower pion masses with the single exponential ansatz. This was a great advantage 
because ansaetze with higher number of exponentials are harder to automate for the error calculation because sometimes the 
order of the masses changes in the fit-results for some bins. We have (almost) always used the single exponential fit for the 



10 



0.065 



0.06 



0.055 



0.05, 



<> 



p 



o 
<:> 



O Fits with from PP 
O Independent fits 



o c c 9 
1> $ l> 



:> 



c;> o 



<> 



10 



12 



FIG. 5: Caa obtained from constrained and unconstrained fits (data slightly shifted horizontally for clarity on the error bars). 



unsmeared correlators both for the pion and the rho. 

Once we determine the pion mass from the PP and the AA correlators, we then respectively use these masses for fitting all the 



correlators to obtain the coefficients, as discussed in Sec. Ill and then extract the various decay constants and the quark masses. 



We have used the same strategy for the smeared propagators too, both for the pion and the rho. 

Fig. [4] shows for K = 0.157 the pion mass determined from single-exponential fits to the PP correlator with local sink and 
smeared source as a function of f m ; n for all the investigated values of the smearing size parameter so = 1,2, . . . , 8. The figure also 
includes the pion mass obtained by single-exponential fits from the local-local PP correlator for this K. The figure shows that 
So = 8 gives the most reliable mass-plateau and we take the mass from the fit with the best confidence level from the plateau at 
So = 8. As discussed below, choice of the optimum so for a given correlator, PP in this case, has to be the same for all smearing 
combinations Is, si and ss. Hence in this case we made sure that for the other two cases, viz., si and ss, the optimum choice of 
so was also 8. Then along with Cf obtained from the correlated 2-parameter single-exponential fits to the Is correlator, we also 
obtained Cf, and C P J at so = 8 by making 1 -parameter single exponential fit with the value of m K put in from the chosen best 
fit to the Is PP correlator. This process was done for each of 200 jackknife bins so as to get the jackknife statistical errors for all 
the derived quantities. The errors shown in Fig. [4] are jackknife errors. 

Once m n is determined from a given correlator (in this Is PP), we can then similarly determine the coefficients for all other 
correlators (in this case, AA, AP and PA) and for all sink-source smearing combinations. Optimum values of the smearing size 
so need not be the same for the other operators. However, it has to be the same for all smearing combinations of the same 
operator (for applicability of Eq. ( 4. 1 1. Once again, for each correlator we looked for the optimum so by comparing 1 -parameter 
single-exponential fits (with m n put in) for all smearing combinations. In Fig. gthe coefficient Cff for the optimum so = 3 for 
the AA correlator is shown as a function of f m j n again at K = 0.157 and displays a nice plateau with accurate data points. The 
figure also compares Cff obtained from correlated 2-parameter single-exponential fit (with no m K predetermined); coefficients 
determined in this way show a rough plateau with a lot less accuracy (larger error bars). This is an important point considering 
that the determination of quark masses and especially the decay constants depend crucially on the accurate determination of the 
coefficients. 



As mentioned already before while discussing the local-local correlators, the pion masses were always determined from only 
the PP and the AA correlators also in the case with smeared correlators. In addition, to be consistent across all K values, we have 
always determined the masses from the local sink - smeared source (Is) combination. This is also true for the rho mass where 
we investigated only the VV correlator. 

For the smeared correlators, we first determine the optimum smearing parameter so for each correlator PP, AA, AP and PA 
and in general the optimum values are different for different operators (see Table |H|). As mentioned already above, for a given 
operator, the value of the optimum so needs to be the same for all smearing combinations Is, si and ss to make use of the formula 
given by Eq. (4.1 1. We have presented all the optimum so values chosen in our analysis for each operator at each K in Table |H| 
There is very little systematics or a visible trend in the values except that usually the PP correlator needed a bigger smearing 
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TABLE II: Optimum values of the smearing size parameter sq for all correlators investigeted 
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FIG. 6: Comparison of Is and si coefficients for all the pion correlators at K = 0.1575. 



size parameter than the AA, the AP or the PA for which one may decipher a rough trend of increasing sq with increasing K. This 
trend is visible also in the VV correlator. 

While the calculation of C.„ is mandatory, in the literature (e.g.. lfL3~ll22l ') C s i (corresponding to smeared sink and local source) 
has often been approximated by C; v and use Cf s /C ss to calculate C//. Since we have performed calculations with both local sink - 
smeared source (Is) and smeared source - local sink (si) combinations, we can check the validity of this approximation . In Fig. 
[6]we plot Cis and C s i at K = 0.1575 for the all four pion correlators as functions of f m ; n . Solid symbols are Is and empty symbols 
are si correlators. As seen in the figure, all coefficients show decent plateau behavior. The data shown are at the respective 
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FIG. 7: Ci[C ss /C s iCi s for all four pion operators as function of K. This ratio should be equal to unity. The blue circles (red squares) represent 
data where the pion mass is determined from the PP (AA) correlator. 



optimum values of the smearing size parameter so for each operator (so — 6 for PP, so = 3 for AA, AP and PA). While for the 
PP and the AA correlators, one should compare their respective Is and si correlators, for the asymmetric correlators one should 
compare Cff with C^f and Cff with Cff. As the figure also shows, it is grossly wrong to approximate Cff by C^f and similarly 
for PA. Actually the difference between Cff and C^ p or between Cff and CS is significant, it is smearing size dependent and 
in this case it is about 50%. However, when one compares the correct quantities, i.e., Cff and C? , , their central values differ by 
about 5%, while the central values of C^f and Cff differ by almost 10%. The difference between Is and si coefficients for the 
case of the AA is the least (about 3%) while for PP it is about 5%. 

In Fig. [7] we show CuC ss /C s iCi s for all four pion correlators as function of K. This ratio should be equal to unity. The circles 
(squares) represent data where the pion mass was determined from the PP (AA) correlator. The straight lines represent average 
values of the respective ratios. They are very close to unity (within a few percent). While the squares (pion mass taken from AA) 
show more fluctuation, the average is almost unity; the circles (mass taken from PP) are more stable, but their average shows 
a few percent bias above unity in all the four correlators. This general character is displayed in all the results of the derived 
quantities we have obtained, viz., data where m n was determined from AA show more fluctuation than the data where m n was 
determined from PP, but show the correct general trend, in particular seem to show less finite size effect (discussed later). 

Before we end the discussion on the details of the implementation of the gaussian smearing, we present two figures to illustrate 
two particular features of the smearing. In Fig. [H]we show an example of oversmearing for the Is PA correlator at K = 0.15675. 
While sq — 1 data look slightly undersmeared (albeit producing a decent plateau), so = 3 completely destroys the plateau, sq = 2 
produces a very stable plateau and has been accepted as the optimum value for this correlator, but the best value could have been 
somewhere between so — 1 and 2 (also please notice the significant difference between the m K values taken from the plateaux at 
Sq = 1 and 2). Our choice of so = 2 for this correlator at this K is also dependent on the proper behavior (i.e. existence of stable 
plateau) of the other two smearing combinations, viz., si and ss of the PA correlator at K = 0.15675. 

Fig. |9]plots ami (computed from Is AA correlator) versus f m ; n at K — 0.1565 for so = 1 to 5. The curious thing to observe in 
the plot is the crossing of the lines approximately at f m j n = 6 indicating that am n is independent of so for the fits using f m ; n = 6. 
One may conclude that the higher state contributions to the correlator are nearly absent at the crossing point and hence the value 
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FIG. 9: Crossing of am n determined from the AA correlator as function of f m j n for different smearing size parameter .so at K = 0.1565. The 
lines are just to guide the eye. 



of am % at the crossing is fairly accurate. This kind of crossing does not take place always, but can be used whenever it does. 
Such a behavior was observed and discussed in the context of static quark potential with gauge field smearing in Ref. j24). 

The analysis for the rho correlator VV is done in a similar way, except that it was less tedious because VV was the only 
correlator we considered. 



VII. RESULTS FOR THE PION: am,, aF n AND am q 



We present in Table III all the results related to the pion. The results are presented in units of the lattice constant a. The 
results are put in two groups depending on whether am K is determined from the PP or the AA correlator, as described in the 



above section (Sec. VI I. In each group, results have been presented for both the local (unsmeared) correlators and the smeared 



correlators. In each of these cases, as explained in Sec III there are three ways each to evaluate aF n and the PCAC quark mass 



14 



0.25 
0.2 

i 

? 0.15 
0.1 
0.05 


0.25 
0.2 

■J 

? 0.15 
0.1 
0.05 




1 


1 1 1 1 1 1 1 

Unsmeared PP 


1 i 


1 i 
























1,1,1,1 








1 


, i 




AA 
AP 
PA 



1 1 


1 1 1 1 1 1 
Smeared PP 


1 


1 1 


i 










Er 






is) 






m 

,i,i,i 








— 

, 1 


1 


□ 

O 

A 


AA 
AP 
PA 



1 1 


1,1,1,1 
Unsmeared AA 


1 


1 1 


i 








& 




® 












□ 


AA 








O 


AP 


, 1 


,i,i,i, 


1 


A 


PA 



_L 



_L 



_L 



_L 



I 1 I 1 I 1 I 1 I 1 I 
Smeared AA 



□ AA 
O AP 

A PA 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.01 0.02 0.03 0.04 0.05 0.06 0.07 



FIG. 10: (antjc) 2 obtained from unsmeared and smeared PP and AA correlators versus different evaluations of am q . 



The errors presented in Table |IH| are statistical errors calculated by the single-omission jackknife method calculated from 200 
jackknife bins. In general, the am K errors increase as K increases. These errors are systematically more for the pion masses 
determined from the AA correlator than from the PP correlator and within each group more for the determination from the 



unsmeared correlators than the smeared correlators. aF n is the most inaccurate quantity presented in Table III In each group, 
both aF K and am q have the least errors when the correlator AP is used; they are the most inaccurate when the correlator PA is 
used. For determination of quantities when pion mass was determined from the unsmeared AA correlator, at times the quality of 
the plateau was poor resulting in inaccurate data; this was especially true for K = 0.15725 and 0.1575. 

plots (am z ) 2 determined from unsmeared and smeared PP and AA correlators versus different evaluations of am q . 



Fig. 
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Except for the data points with the lowest pion and quark masses, the rest appears to be roughly consistent with the lowest order 
(LO) chiral perturbation theory (%PT). As already mentioned, out of the four plots in Fig. 10 the bottom left plot (pion mass 
from unsmeared AA) looks the worst. In all cases, the PCAC quark mass am q determined from the AP and the PP correlator 
(shown in plots as AP: open circles) has the least error. 

Any departure from LO %PT shows up in the behavior of the ratio (am K ) 2 / (am q ), called the chiral ratio below, as a function 
am q . This is plotted for the data obtained only from the smeared correlators in Fig. 



1 1 Firstly we need to separate any effect of 



finite size effect on the masses before discussing the chiral behavior. We are quite confident that for the values of K from 0. 156 to 
0.1575, there is negligible finite size effect on these masses, because data, consistent with our numbers, are available in previous 
literature at some of these K values at larger volumes (see below for discussion and references), For the masses at K = 0.15775, 
our guess is that for different operators there is a difference in finite size effect. From both the Figs. 10 and 11 it seems that 



there is a little bit of finite size effect for the pion mass at this K, but more with pion mass determined from PP than AA, because 



in the upper plot of Fig. 1 1 there is an upward trend at this point. This is not to be taken as an effect from next-to-leading order 
(NLO) %PT because with standard values of low energy effective (LE) constants of %PT obtained from other studies (both on 
the lattice and otherwise), there should not be an upward trend of the ratio (as seen in the upper plot here) at this pion mass. 
Difference in the finite size effect on the pion mass calculated from different operators is an interesting observation and to the 
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TABLE III: Results for am n , aF % and PCAC quark mass am q . 
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FIG. 11: The ratio (am,) 2 /( am„) obtained from smeared PP and AA correlators versus different evaluations of am a 



best of our knowledge has not been discussed in the literature before. It is fair to say that this is not an unexpected behavior. 
However, what is also interesting is that the generally accepted noisier operator AA produces less finite size effect. 

If our data see only the LO %P1 ', then the ratio should be a perfect straight line parallel to the am q axis. If we believe our 
error bars, there is a significant downward trend of the ratio from the larger quark masses to the smaller quark masses (more 



pronounced with the AA and PA quark masses). For the lower plot of Fig. 1 1 this trend continues and intensifies to the second 



last point from the left (at K = 0. 15775) which can only be an effect from higher order terms of %PT . 

The left-most point (at K = 0.158) definitely has a significant finite size effect as we can compare our data with evaluations on 
larger volumes (see below). Even at this point, the pion mass determined from the smeared AA correlator has smaller finite size 
effect as apparent from the figure. We also like to point out that the finite size effects, its dependence on the pion operators etc 
have nothing to do with smearing or the smearing size [ 19 1. The same effects are also visible from similar data from unsmeared 
operators (although not plotted in Fig. [TT] but available in Table |HI| and Fig. [T0| . 

Fig. 12 plots aF K versus am q in way similar to Fig. 10 The bottom left plot where the pion mass is determined from unsmeared 
AA correlator has the most inaccuracies. The main difference with the am\ plot is that the finite size effect at K = 0.15775 and 
0.158 is now quite severe. This is consistent with %PT which predicts four times larger (and opposite in sign) finite size effect 
for aF K as compared to arn^ [25 1. The behavior of aF K with am q using the AA operator (squares) show the most continuous 
behavior, even though these have larger error-bars compared to the am q data computed from the AP correlator. 
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FIG. 12: aFn obtained from unsmeared and smeared PP and AA correlators versus different evaluations of am q . 



A. Chiral behavior 



We have plotted in Figs. 13 and 14 the whole content of the bottom figures of Fig. 1 1 and[T2]respectively (data from smeared 
AA correlator only, represented by open symbols) along with (unimproved) Nf = 2 Wilson data at — 5.6 from Ref. Il22l 
(16 3 32 and 24 3 40), Ref. |26| (16 3 48) and Ref. £Q (24 3 32) (data from others represented by different filled symbols). These 
figures also show NLO %PT El plots 
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with aF — 0.035, aB — 1.5. With these values of aF and aB as input we actually could fit our data (for amj^, represented by 
squares) for both the chiral ratio and aF n . For the chiral ratio, the fits include points from K = 0.15675,0.157,0.15725 and 
0.15775. The point at K = 0.1575 is excluded from the fit because of systematic fluctuations. The fit to the aF K data is also 
done with the AP and PA quark masses in a similar range of data. The values of aA 3 and aA^ that come out of the fits are: 
aA 3 = 0.374(7), aA 4 = 0.813(11). 

Let us state categorically that the chiral NLO plots are not to be taken as serious fits, they are shown more to understand our 
data and where the data of other works (with same parameters), especially on larger volumes, are in relation to our data and the 
chiral NLO plots. Firstly, the NLO fits for both the chiral ratio and the aF K are done with the same input values of aF and aB. 
Secondly, the fit to the chiral ratio is done with amj^ (denoted by open squares) and nearly goes through all the points of the 
fit, Our other data points with am q p (circles, with their small error bars) and with am PA (diamonds, with large error bars) are 
also close to this NLO curve. Thirdly, the data points on larger lattices (at K = 0.158 by J4] 122), at K = 0.1575 by |4| lie very 
close to both the NLO fits (Figs 



0.15825 has a large error 

bar on the chiral ratio plot, however, on the aF K plot looks very close to the NLO fit. The errors for the 16 3 data from l22"ll26l 



13 and[14|i. The data with the smallest quark mass from [4| at K - 
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FIG. 13: The chiral ratio (anix) 2 / am q obtained from smeared AA correlator versus different evaluations of am q . Figure includes data from 
other work at same or larger volume and a suggestive NLO chiral plot. 



are calculated by us from their quoted pion and quark mass data by quadrature (Ref. [26 1 has no F K data either). Some of these 
16 3 data look a bit erratic, especially from ll22l . Finally, probably because of large finite size effect on dF % , the numerical data 
computed on relatively smaller lattices (including ours) tend to have a larger downward bend than is consistent with acceptable 
values for the low energy constants. However, our fits with aF and aB as inputs, done from K = 0.1565 to 0.1575, for all the 
quark masses am^, am q p and amj 4 give a reasonable fit with the evaluations upto K = 0.15825 on larger lattices staying very 
close to this explorative NLO fit. 

At this point if we look back at the Fig. 10 and take a closer look, we discover that this figure is not consistent with LO %PT, 
because data at larger masses do not go though the origin in a straight line. The data in all the plots of Fig. 
small bend. 



10 actually have a 



B. Comparison of VWI and AWI masses 

In all the Tables and Figures above, we have presented results for the quark mass from the PCAC relation, am^ WI , which 
is devoid of &{a) effects because the PCAC relation guarantees that the square of the pion mass vanishes as the quark mass 
approaches zero. %PT can also be made consistent with this, once the quark masses used there are the PCAC masses |28|. 

Traditionally, one also defines the VWI (vector Ward Identity) quark mass as 

where K — > K c = 1/8 = 0.125 in the chiral limit for the free theory (<7 X/J — 1). 

At the tree level, this definition takes care of the &{a) additive quark mass. However, at the quantum level, ni^ WI should be 
devoid of 6(a) effects only approximately. The effect should be worse at smaller values of j5 with a larger lattice constant a. 
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FIG. 14: aFx obtained from smeared AA correlator versus different evaluations of am q . Figure includes data from other work at same or larger 
volume and a suggestive NLO chiral plot 
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FIG. 15: (amx) 2 obtained from unsmeared and smeared PP (left) and AA (right) correlators versus 



The behaviour of (am„) 2 as function of ^ is shown in Fig. 15 for pion masses determined from unsmeared and smeared 



PP and AA correlators. Although the data from unsmeared AA correlators show more fluctuations, straight lines fits excluding 
K = 0.158 (because of sizable finite size effects) are possible in all cases giving a determination of K c . The values of K c obtained 
are as follows: (i) 0.15858(3) (PP-smeared), (ii) 0.15862(4) (PP-unsmeared), (iii) 0.15854(4) (AA-smeared), and (iv) 0.15855(6) 
(AA-unsmeared). For the same lattice volum e and j3, SESAME and TjL Collaboration |29| obtained K c = 0.158493(16). 
The values of am^ WI calculated using Eq. (7.3 1 are presented in Table IV 
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K 


ami from PP 


am n from AA 


with smearing 


without smearing 


with smearing 


without smearing 


0.156 


0.05215(55) 


0.05294(75) 


0.05135(75) 


0.05155(126) 


0.1565 


0.04191(55) 


0.04270(75) 


0.04111(75) 


0.04131(126) 


0.15675 


0.03681(55) 


0.03761(75) 


0.03601(75) 


0.03621(126) 


0.157 


0.03173(55) 


0.03253(75) 


0.03094(75) 


0.03113(126) 


0.15725 


0.02667(55) 


0.02746(75) 


0.02587(75) 


0.02607(126) 


0.1575 


0.02162(55) 


0.02242(75) 


0.02083(75) 


0.02102(126) 


0.15775 


0.01659(55) 


0.01738(75) 


0.01579(75) 


0.01599(126) 


0.158 


0.01157(55) 


0.01237(75) 


0.01078(75) 


0.01098(126) 



TABLE IV: am vwl in dependence of K. 
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FIG. 16: Deviation of of m^ WI from n£ WI (some of the data slightly shifted horizontally for clarity). 



We note that has no 6(d) scaling violation whereas m q has some remaining 0{d) scaling violation. Fig. 16 plots the 



deviation of m% WI from mf™ 1 as a function of i. We notice small fluctuations (a few percent) around a normalization factor of 
approximately 0.78 (i.e., m^ WI w 0.78/nY ). Even at K = 0.158 where the fluctuation is maximum shows a deviation of about 
6% from the normalization. If these deviations are any signature, scaling violations may be small in our case, This is to be 
expected at our lattice scale a « 0.08fm (aT 1 ~ 2.45GeV) 0. 



VIII. RESULTS FOR THE RHO: am p , aF p AND DIFFERENT RATIOS 

Table |Vj presents our results for the rho mass and the decay constant in lattice units for all the K values of our simulation. 
The errors are again the single-omission jackknife errors computed from 200 jackknife bins. In general the errors are more for 
larger K and for data from unsmeared correlators (especially for am p ). With the range of pion masses reached in our simulations, 
2m K > m p was always satisfied and there was no complication in the investigation of the rho correlator. 

The data for am p is plotted in Fig. lT7las a function of am^ p . For clarity of the figure, other quark masses are not shown in this 
figure (even in this case, there are fouraifferent PCAC quark masses depending on whether the pion mass was determined from 
unsmeared or smeared PP or AA correlators). The straight line fits shown in the figure exclude points for the two extreme K and 
for K = 0. 15725 (appearing to have a systematic error) and lead to a value close to the physical rho mass on chiral extrapolation. 
As expected, the lowest rho mass at K = 0.158 has some finite size effect. 



Fig. 18 shows three ratios F n /F p , 1 / f p = F p /m p and F n /m p obtained from smeared correlators versus PCAC quark masses. 
The left (right) figure shows the dependence of these ratios on AA, AP and PA quark masses when the pion mass was determined 
from the smeared PP (AA) correlator. The drop at small quark masses in the ratios F n /F p , and F n /m p are mainly attributable to 
the most significant finite size effect of F K . The quantity in the middle 1 / f p = F p jm p is actually a ratio of the other two ratios 
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K 


am p 


aF p 


with smearing 


without smearing 


with smearing 


without smearing 


0.156 


0.5365(23) 


0.5409(32) 


0.1577(17) 


0.1631(24) 


0.1565 


0.4935(25) 


0.4890(41) 


0.1460(17) 


0.1443(27) 


0.15675 


0.4730(27) 


0.4765(43) 


0.1427(16) 


0.1460(28) 


0.157 


0.4516(28) 


0.4510(50) 


0.1369(13) 


0.1385(30) 


0.15725 


0.4170(45) 


0.4202(56) 


0.1266(31) 


0.1300(30) 


0.1575 


0.4069(34) 


0.4098(70) 


0.1290(25) 


0.1323(40) 


0.15775 


0.3811(70) 


0.3847(82) 


0.1175(47) 


0.1221(44) 


0.158 


0.3675(71) 


0.3697(122) 


0.1103(99) 


0.1130(62) 



TABLE V: am and aF from smeared and unsmeared VV correlators for all the values of K 
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FIG. 17: amp obtained from unsmeared and smeared PP and AA correlators versus am^ p . 



and with the elimination of F n cancels out whatever finite size effects of F p and m p . 

With phenomenological inputs of m p = Hl.lMeV , F K = 92AMeV and y/2F p = 2\6MeV, Fig. 18 also enables approximate 
estimations of the renormalization constants Zy and Z& associated with F p and F n respectively. Fit to all points of the ratio 
l//p = F p /m p leads to Zy = 0.659(2) and a similar fit to points corresponding to K =0.1565 to 0.15725 for the ratio F K /m p 
leads to Za = 0.76(1). These are remarkably close to the values quoted in Il3l l26l . 

IX. SUMMARY AND CONCLUSIONS 

Due to recent theoretical, algorithmic and technological advances, Lattice QCD with Wilson fermions has entered a very active 
phase. We have planned on a detailed lattice QCD investigation with Wilson and Wilson-type fermions. In this work we have 
taken the first step towards addressing the issues that are encountered probing the chiral regime of lattice QCD systematically. 
Here, our emphasis is to have fully dynamical simulations at an extensive set of K values and study accurate determinations of the 
masses and decay constants in as many ways as possible. We employ standard (unimproved) Wilson gauge and fermion actions 
with 2 fully dynamical light quark flavors on 16 3 32 lattice, j3 = 6/g 2 = 5.6 at 8 values of bare quark masses (corresponding to 
fermionic hopping parameter K= 0.156, 0.1565, 0.15675, 0.157, 0.15725, 0.1575, 0.15775, 0.158). The lattice scale reached 
is respectable and a ~ 0.08 fm (cr 1 ~ 2.45 GeV) @. 

In order to accurately extract the mass gap from the asymptotic behavior of Euclidean time correlation functions, various 
smearing techniques applied to the hadron operator have been employed in the literature. Among them, gaussian smearing 
which can be motivated from conceptually simple considerations can be implemented in numerical simulations in a simple and 
straight forward manner. To assess the strengths and weaknesses of the gaussian smearing, we have carried out an extensive 
study that employs local-smear, smear-local and smear-smear sink-source combinations at many values of the smearing size 
parameter sq. 
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FIG. 18: F x /Fp, 1/ fp = Fp /trip and F K /nip obtained from smeared correlators versus different evaluations of am q 



At each value of the smearing size, to extract meson observables, a variety of correlation functions are measured. Using the 
pseudoscalar and the fourth component of the axial vector densities denoted by P and A respectively, for the pion, PP, AA, AP 
and PA correlators are used since both P and A carry quantum numbers of pion. We have extracted mass gaps and amplitudes 
(coefficients of the correlators) for the smearing radius so ranging from 1 to 8 in steps of 1 for all these corelation functions for 
all the three smearing combinations. At each K, we find the optimum smearing size sq for each correlator type; the optimum so 
in general is different for different correlators types and also varies with K. In the case of pion, we do not detect any systematics 



in the optimum choice of sq as a function of K as is seen from Table III except that the optimum smearing size is smaller for AA, 
AP and PA than for PP. One should also be cautious about oversmearing. In the case of rho we find that sq gradually increases 
as K increases. 

We have also carried out extensive calculations with local operators in order to quantify the systematic effects in the determi- 
nation of masses and decay constants. 

We have obtained the pion masses from PP and also from AA correlators. These masses are then used to determine the 
coefficients for all the correlators. In each case we obtained pion decay constant and quark masses using AA, AP and PA 
correlators. The values obtained using AP (PA) correlator usually has the lowest (highest) errors. 

The general conclusion regarding the use of the different correlator types PP, AA, AP and PA appears to be the following: 
whenever the axial vector operator is at the source (i.e., AA or PA), the data (be it the pion mass, decay constant or the quark 
mass) seems to have more noise and statistical errors. However, despite that, these data seem to have a better general trend 
overall (be it consistency check for smearing (Fig. 17), comparison with the VWI quark mass (Fig. 16 1 or chiral trends (Figs. 13 
and 14 1, they also seem to have significantly less finite volume effects. 

The quantities derived from correlators with the psedoscalar operator (i.e., PP and AP) at the source seem to have the least 
statistical errors, however, they seem to carry some systematic errors and definitely show signs of stronger finite size effects at 
larger K values. 

Although we have performed a lattice QCD investigation with fully dynamical quarks at a host of values of K, we do not 
have nearly enough data points or small enough quark masses to have a reliable chiral extrapolation. Let us make it very clear 
that although it is one of our ultimate aims, in this paper it was not our intention either. Given our emphasis on other aspects, 
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e.g., study of smearing in a very detailed manner or comparison of all pion correlators, we still can make a few interesting 
observations regarding chiral trend of our data. In both the behaviors of aF K and the chiral ratio {am n ) 2 / (am q ) versus am q , we 
see indication of departure from LO %PT . The chiral ratio has a slope downwards in the range of am q = 0.05 - 0.025 and our data 
is more consistent with NLO %PT than LO. Data from other works (with same parameters as ours) done at larger volumes are 
also consistent with this chiral trend including points at smaller quark masses upto am q < 0.015. The chiral ratios at K = 0.158 
with 16 3 32 lattices (from all collaborations including ours) show significant finite size effect. 

The data for aF K is also consistent with NLO %PT in the range am q = 0.05 - 0.03. The faster decrease with respect to am q 
at K = 0.15775 and 0.158 is clearly due to, and in accordance with %PT finite size effects. This is again corroborated with the 
simulations done at the same and larger volumes for the same K range. 

We have studied the deviation of am v q WI from am^ WI as a function of K. We note that has no &{a) scaling violation 

whereas am™ 1 has some remaining &(a) scaling violations. We have noticed small fluctuations (a few percent) around a 
normalization factor of approximately 0.78. If these deviations are any signature, scaling violations may be small in our case. 

The dimensionless ratio -f- is remarkably independent of am q in the range 0.07 - 0.015 and we can easily extract the constant 
value of this ratio 0.301. Comparison with the physical value of this ratio (0.198) yields an estimate of Zy « 0.659. The 
dimensionless ratio is approximately independent of am q in the range 0.055 - 0.03. The deviation from constancy for the 
smallest two values of am q is due to finite size effects. The constant value extracted for the ratio is 0.159. Comparison with the 
physical value of this ratio (0.120) yields an estimate of Za ~ 0.76. It is encouraging that the extracted values for both these 
renormalization constants (which are independent of the lattice scale) are in agreement with those extracted non-perturbatively 
on the lattice in Ref. ||3l l26l . 

We find our initial results interesting and very encouraging. Simulations are underway utilizing significantly larger lattices 
and improved algorithms in order to further test the accessibility of chiral region to the naive Wilson action (both gauge and 
fermion sectors). It will be interesting to compare our results obtained with unimproved Wilson fermion and gauge action to that 
with improved actions, for example, Refs. |[30"ll3T1[32l . 
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APPENDIX A: FOURIER TRANSFORM TRICK FOR SINK SMEARING 

In this appendix, we follow the discussion in Ref. [ 16 1. Generically we need to calculate 

£(^ slnk (x,f)^ ource (x ,r )) . (Al) 

x 

Case I (local source, local sink) 
We have 

£(V7(x,f)rv/(x,f)v7(xo,fo)rv/(x ,fo)) = £(Trace (S(x,f;x ,fo)rS(x ,f ;x,f)r)) (A2) 

X X 

where T denotes any of the sixteen Dirac matrices. The fermion propagator S(x,y) is obtained by solving the equation 

MS =8 (A3) 

where M is the Dirac- Wilson operator. 

Writing the coordinates explicitly, we have, 

£ M(zi,T;xi,fi) S(xi,fi;yi,0) = <5 Zl , yi 5 T , . (A4) 



For the calculation of the meson propagator, Eq. (A2i, we need to calculate only S(x,f;xo,fo) since S(xo,fo;x,f) is trivially 
calculated from 

S(x ,fo;x,f) = y 5 S\x,t;x ,t ) Ys ■ (A5) 
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Case II (local source, smeared sink) 



^sink 

XiXt 



£ 0(X1-X) 0(X 2 -X) V(xi,f)r l^(X2,f) 
x l x 2 

^source = . (A6) 

Now we have, 

£ £ 0(Xi -x) <j>(x 2 -x) (Trace 5(x 2 ,f ;x ,f ) r7 5 5 t (x 1 ,f;x 0! fo)75 . (A7) 

X X[X2 

Going to momentum space, 

0(xi -Xo) - 1/iV £ e -'kr("i-x ) ^ kl ^ 

ki 

5(x 2 ,f;x ,fo) - £ «-*3-(»z-»o) S(k 3 ,0 , etc., (A8) 

k 3 

we get, 

(1/^V) 4 £ 0(k) 0(k) (Trace S(k,f) T y 5 §\k,t) Y 5 T) . (A9) 

k 

Need to calculate only 5(k,f) and 0(k). 
Case III (smeared source, local sink) 

^source = £ - *o) <M*2 ~ X ) y7(*l ,'o)r V/(x 2 ,f ) 

X[X 2 

^wk = jjf(x,t)r \j/(x,t) . (A10) 



Now we have, 



X X^X2 



where 



£ £ 0(xi-xo) 0(x 2 -x o ) (Trace5(x,f;x 1 ,f )r755 t (x,f;x 2 ,fo)7 5 r) . (All) 

= £ (Trace^(x,f;x ,fo)r75^ t (x,f;x ,f )75 T) (A12) 

X 

y(x,t;xo,t ) = £ ^(x 1 -x )S(x,t;x 1 ,t ) . (A13) 



Starting from Eq. ( |A4| , 



£M(y,f,;x,f)5(x,f;x 1 ,0) 0( Xl -x o ) = £5(y-xi) 0(xj -xo) . (A14) 

X] X, 

i.e., 

M(y,f v ;x,f) J^(v;x ,0) = </)(y-x ) S(* y ) . (A15) 

Once 5^ is calculated, is calculated trivially. 

Note that in this case there is no need to Fourier Transform. 

Case IV (smeared source, smeared sink) 

In this case again it is advantageous to perform Fourier Transform and we get 

£(^ sink (x,f)^ ource (x ,fo)> - £ m 0(k)(^(k, f ) r 75 ,y\Kt) r5 r> . (Ai6) 

X k 
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